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Abstract 



CD ■ Chaotic dynamics can be effectively studied by continuation from 

, an anti-integrable limit. We use this limit to assign global symbols to 

' orbits and use continuation from the limit to study their bifurcations. 

^\ , We find a bound on the parameter range for which the Henon map 

exhibits a complete binary horseshoe as well as a subshift of finite 
type. We classify homoclinic bifurcations, and study those for the area 
'"^ ' preserving case in detail. Simple forcing relations between homoclinic 

^ , orbits are established. We show that a symmetry of the map gives rise 

' to constraints on certain sequences of homoclinic bifurcations. Our 

, numerical studies also identify the bifurcations that bound intervals 

on which the topological entropy is apparently constant. 
AMS classification sclieme numbers: 58F05, 58F03, 58C15 

• I— I 

X. 

^ . 1 Introduction 

The problem of determining the sequence of bifurcations that result in the 
creation of a Smale horseshoe is an interesting one ||l|, ^. In this paper 
we use a continuation technique based on an "anti-integrable" (AI) limit 

to study some of these bifurcations from the opposite side, that is, as 
bifurcations that destroy the horseshoe. 

As a simple example, we study the family of Henon maps js], |6| 

x'\_fy-k + x^\ , . 

y' )-\ -bx )■ 

Apart from the fact that the Henon maps are the simplest, non-trivial maps 
of the plane, they are of more general interest as well, since vector fields 
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in the neighborhood of certain co dimension- two homochnic bifurcations can 
be reduced to Henon-hke maps 0, 

As we recall in §|2|, the AI limit for this map is essentially k ^ oo. In 
order to represent this with finite parameters, we need only rescale the map, 
letting 

z = ex , where e = —= . 

Vk 

As was shown by Devaney and Nitecki |^ , the Henon map has a hyper- 
bolic horseshoe when 

^>(l + l&l)^^^. (2) 



The Henon map has at most 2" ||10[] periodic points of period n, and when 
the map has a hyperbolic horseshoe, all these orbits exist and are easily 
identified by their symbolic labels. 

We showed earlier |11] that a contraction mapping argument implies 



there is a one-to-one correspondence between orbits in the AI limit and 
bounded orbits of the Henon map in precisely the range Eq. (^). In §|5| we 
show that if we consider a particular subset of orbits, this bound can be 
increased. Moreover, in §p, we present results of numerical investigations 
for all 6 that give what we believe are optimal bounds. 

In general, the existence of an anti-integrable limit leads to a natural 
symbolic characterization of orbits — for the Henon map this is the same as 
the horseshoe coding. We use this coding, and as discussed in §|3|, a predictor- 
corrector continuation method |12|, to give each orbit a global code. That is. 



we label an orbit with the AI code, and use this designation for the family 
of orbits until it collides with a family with a different code. For the Henon 
map, this gives a map from the bounded orbits of the map to sequences of 
symbols s S 2^ modulo cyclic permutations, providing only that every orbit 
can be smoothly connected to the AI limit. This is the working hypothesis 
for our numerical method, even though we know that it is probably not true 
in general. It is certainly valid when the hyperbolic horseshoe exists. We 
give an example of a periodic orbit not smoothly connected to the AI limit in 
the dissipative case. This illustrates a general anti-monotonicity result p^ , 
stating that when b ^ —1, 0, 1 the map generically creates but also destroys 
orbits when k is increased, so that the topological entropy is not necessarily 
monotone. We extend this anti-monotonicity result to the area preserving 
case by a quite different argument concerning the vanishing of twist in the 



neighborhood of the period tripling bifurcations in a separate paper |14]. 
Even though the area preserving case exhibits anti-monotonicity, we still 
conjecture that there are no isolated bubbles in the bifurcation diagram, i.e. 
that every orbit is continuously connected to the AI limit. 

Our global code contrasts with other methods for constructing symbolic 
dynamics for maps, which rely on some attempt to obtain a generating parti- 
tion ||T5| , |2|, IC, 17, 18]. These methods rely on somewhat ad hoc techniques 
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for constructing the partition, especially when there exist elliptic orbits. 
Our method gives a natural partition that is smoothly connected to the 
horseshoe, though it does rely on our working hypothesis 

In our computations of the Henon map, we observe that the horseshoe 
destroying bifurcation appears to be a homoclinic saddle-node bifurcation 
when the map is orientation preserving, and a heteroclinic saddle-node when 
it is not. In §0 we study in detail the homoclinic orbits of the area preserving 
Henon map, and show how the AI code directly gives other properties of the 
orbits, such as their "type," "transition time," and "Poincare signature." 

For an area preserving map, the destruction of a horseshoe by a homo- 
clinic bifurcation gives rise generically to elliptic orbits. Specifically, if / is 
a area preserving diffeomorphism with a homoclinic tangency at x then 
for any neighborhood U of x, there is an area preserving diffeomorphism 
C^-close to / that has an elliptic periodic point in U [19|. 

Much more is known about the behavior of area-contracting maps near 
a homoclinic tangency. Gavrilov and Silnikov proved that if a map 
has a quadratic homoclinic tangency at a parameter k* then there exists a 
sequence of parameter values A;„ — > k* such that at kn there is a saddle-node 
bifurcation creating orbits of period n [ pO| , 21 1; because one of the created 
orbits is a sink, this is called a cascade of sinks. Robinson extended these 
results to the real analytic case where the intersection is created degenerately 
[p2| . In our computations we will find a similar cascade of saddle- node 
bifurcations for the area preserving Henon map — this gives a sequence of 
elliptic orbits that limit on the homoclinic bifurcation. Thus the destruction 
of the horseshoe is associated with the creation of the first stable "island." 

The ordering on the invariant manifolds poses severe restrictions on the 
possible bifurcations. In §P we use it to prove which homoclinic bifurcation 
of the hyperbolic fixed point is the first one. We observe that the forc- 
ing relations between homoclinic orbits up to type 6 is essentially like the 
unimodal ordering of one dimensional maps. Generically a homoclinic bi- 
furcation corresponds to a quadratic tangency of the stable and unstable 
manifolds — a "homoclinic saddle-node bifurcation." There are two more 
generic bifurcations in maps with a symmetry: a homoclinic pitchfork when 
the manifolds exhibit a cubic tangency, and a simultaneous pair of asym- 
metric saddle node bifurcations. In §^we show that a symmetric homoclinic 
bifurcation forbids certain other bifurcations to occur after it, leading to a 
natural mechanism to create homoclinic pitchfork bifurcations or asymmet- 
ric saddle node pairs. We observe all three of these bifurcations for the area 
preserving Henon map, which has a time-reversal symmetry. 

Davis, MacKay, and Sannami conjectured that there are intervals of k 
below the horseshoe for which the Henon map is a hyperbolic Markov shift. 
They also identified the Markov partitions for these cases. Their conjecture 
was based on computing all the periodic orbits up to a period 20 using the 
technique of Biham and Wenzel [23, 24|. In §^, we confirm their compu- 
tations with our continuation technique and extend them to period 24 — an 
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order of magnitude more orbits. Moreover, we identify the bifurcations re- 
sponsible for the creation and destruction of these apparently hyperbolic 
intervals; as befits with the theme of this paper, they are homoclinic bifur- 
cations. 



2 Anti-Integrable Limit 

Dynamics in discrete time can be represented by a relation, F(x, x') = 
where x and x' are points in some manifold. Normally, we can explicitly 
solve for x' = f{x), giving a map on the manifold, with orbits defined 
by sequences xt = f{xt-i). Suppose, however, that F depends upon a 
parameter e, in such a way that this is not always possible; for example, 
F[x,x') = eG{x,x') + H{x). In this case the implicit equation F = 
can no longer be solved for x' when e = 0; instead "orbits" correspond to 
arbitrary sequences of points, xt that are zeros of H — the dynamics are 
not deterministic. In this case we say that e = corresponds to an anti- 
integrahle (AI) limit of the map / [^5|. If the derivative of H is nonsingular, 
then a straightforward implicit function argument can be used to show that 
the AI orbits can be continued for e 7^ to orbits of the map / [Q, 26 1. An 
AI limit with this property is called nondegenerate. 

For example, consider the Henon map Eq. (||). Denoting points on an 
orbit by a sequence xt, t £ Z, we can rewrite Eq. (||) as a second order 
difference equation 

xt+i + bxt-i + k - xf = . 

Introducing the scaled coordinate z = ex and choosing e = k^^^'^ gives an 
implicit map in the variable z with parameter e 

e(zt+i + 6zt_i) + 1-^2 = 0. (3) 

With this choice of e, we can study only the range < A; < 00; however, 
one could redefine e to shift this range.|^ A period n orbit of the Henon map 
is given by a sequence zo,zi,. . . ,Zn-i that satisfies Eq. @, together with 
the condition that Zf+n = zt- The corresponding family of periodic orbits is 
denoted by z(e). 

At the AI limit, the map Eq. (js]) reduces to 

z^ = l. 

Thus "orbits" in this limit are arbitrary sequences of ±1, which we abbre- 
viate with + and — . Let S denote the space of such sequences 

S = {-,+}^ = {s : St G {-,+}, tGZ}. (4) 



^ For example, choosing e = (fc + 5)"^^^, maps positive values of e to the range ~5 < 
k < 00. Our numerical routines typically use 5 = 1 so that we can cover the entire 
parameter range where there are bounded orbits. In the text we always use 5 — 0. 
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For ease of notation we denote the corresponding sequence of {+1, —1} G 
by the same symbol s. Hence every sequence s G S is an orbit s G at 
the anti-integrable hmit, and each of these can be continued to an orbit of 
the Henon map for smah enough e [27, Previously we gave an explicit 
upper bound on e for the existence of orbits for every symbol sequence |11]: 

Theorem 1. For every symbol sequence s G S there exists a unique orbit 
z(e) of the Henon map Eq. such that z(0) = s providing 



e|(l + |6|) <2v/l-2/\/5 0.649839. (5) 



The basic idea of the proof of this theorem is as follows |11]. Let Bm be 
the closed ball of radius M around the point s G S, 

5A/(s) = {z:||z-s||oo<Af}. (6) 

For each symbol sequence s G S and small enough M, define a map T : 
Bm Bm by 



Tj(z) = Si'\/l + e{zi+i + bzi^i) , (7) 

then the corresponding Henon map orbit z(e) is a fixed point of T. The 
conclusion of Theorem |l] follows from finding the maximum value of e for 
which there is an n such that T" is a contraction mapping (i.e., T" : Bm 
Bm and ||DT"|| < 1). 

The fact that T is a contraction implies that there are no bifurcations 
in the range Eq. (||). This is the statement of 

Corollary 2. There are no bifurcations in the Henon map when e and b are 
in the range given in Theorem 

Proof: Denote the system of equations (^) by H(z, e) = 0. This infinite 
continuation problem has a unique solution z{e) if the inverse of -D^H is 
bounded. The t-th component of T is related to the t-th component of H 

by 

Ht = - zl 

Differentiating this at the fixed point that exists according to Theorem || 
gives 

L>3 = diag(2zt)(DT-id). 

The operator has bounded inverse because zt is bounded away from 

zero by Eq. (^) and the inverse of DT — id exists because ||-DT|| < 1. □ 



This result can be extended to imply that the invariant set is uniformly 
hyperbolic for the case 6 = 1 when the operator D^H is symmetric [p8| . 
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In §11 we will use numerical continuation to estimate the parameters at 
which the first bifurcation occurs, giving an improvement in this bound, 
albeit a numerical one. 

It is interesting that Devaney and Nitecki Q obtained precisely the same 
bound, Eq. (^), for the parameter domain in which the non- wandering set 
of the Henon map is a hyperbolic horseshoe. Nevertheless the AI contin- 
uation argument has has two advantages over the geometrical arguments 
of Devaney and Nitecki. First, it easily generalizes to higher dimensions 
allowing one to compute parameter bounds for the existence of horseshoes 
in higher dimensional maps |29]. Second, it allows us to easily bound the 



parameter range for which certain subsets of orbits exist (i.e. the parameter 
range where the map is conjugate to a subshift of finite type). We present 
such a bound for a subshift of finite type in §|5|. 

3 Numerical Method 

In this section we formulate the problem of following Henon map orbits away 
from the anti-integrable limit as a classical continuation problem pO| . 

A period n orbit family of the Henon map with coordinates given by z(e) 
is a zero of the function G : M" x M — > whose t^^ component is given by 
the left hand side of Eq. @. The zeros of G are generically smooth curves 
in M" X M defined by the continuation problem 

G(z,e) = 
z(0) = s ■ 

Practically the continuation is always started at the AI limit. The curve 
might either extend to some maximal emax and return the the limit or just 
continue indefinitely. Since for the Henon map there are no orbits for k < —1 
we expect that all the curves do return to the AI limit. 

This is a standard continuation problem, which we solve using a predictor- 
corrector method with a linear tangent predictor and a Newton's method 
corrector. For numerical linear algebra we use the Meschach library |31]. 



The algorithm incorporates an adaptive step size control with bisection 
backtracking if the corrector fails to converge. The algorithm terminates 
when a user-specified value of e is reached or, when the tangent direction 
is not uniquely defined. The process of continuing a sequence of orbits is 
trivially parallelizable since the operations performed on each orbit are com- 
pletely independent of each other. We use a "divide and conquer" strategy 
to spread the total computational effort across several different machines 
running simultaneously. This is especially advantageous when the number 
of orbits continued reaches into the millions. 

Continuation methods are based on the assumption that the orbits of 
interest are actually connected to the limit at which the continuation starts. 
Since the Henon map does not have an integrable limit, the natural starting 
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point would be to continue all the orbits that bifurcate off the fixed points. 
But this would only yield a small fraction of all orbits: many of them are 
born in saddle-node bifurcations that are not connected to any other or- 
bit. Therefore the AI limit is a much better limit from which to continue. 
However, the same general restriction applies, i.e. only orbits that are con- 
nected to the limit (or their parents, grandparents etc.) can be found. We 
formulate this central hypothesis as the "no-bubble-conjecture": 

Conjecture 1 (No Bubbles). For the area preserving Henon map every 
orbit is (at least) continuously connected to the anti-integrable limit. 

We could only hope for continuous connection in M" x M because the 
branches corresponding to the children in a bifurcation are in general not 
smoothly connected to the parent. Our numerics currently does not perform 
any branch-switching from parents to children. Therefore, in practice, we 
are actually using the working hypothesis: "Every orbit of the Henon map 
can be smoothly connected to the anti-integrable limit." Unfortunately our 
working hypothesis is not true in general. For b ^ —1,0,1 it has been 



shown [13 1 that periodic orbits are both created and destroyed when the map 



parameter k is increased (the authors of |]1^] call this "antimonotonicity" ^ 



Moreover, in [14| we show that this conclusion holds for the cases 6 = ±1 
as well. Consequently the topological entropy is not necessarily a monotone 
increasing function, as it is for the logistic map, 6 = 0. In the following we 
will elucidate the relation between antimonotonicity, our working hypothesis 
and the no-bubble conjecture. 

Even though we have antimonotonicity whenever b ^ 0, this does not 
readily imply that our working hypothesis is false. In particular the smallest 
period orbit that is antimonotonic when 6 = 1 is of period 10, and it is still 
smoothly connected with the AI limit [U]. However, orbits that bifurcate 
from this one may violate the working hypothesis 

The worst possible case from the point of view of continuation is an 
orbit that neither smoothly nor continuously connects to the limit, i.e. an 
"isolated bubble." Note that in the area preserving case, one orbit of this 
type implies an infinite number of them because it must be born in a saddle 
node bifurcation and the stable orbit of the created pair generically passes 
through an infinite number of rational winding numbers. In order to con- 
stitute a violation of our conjecture, none of these orbits would be allowed 
to reach the AI limit; otherwise the original orbit would be continuously 
connected. 



^ This orbit is created in a 3/10 rotational bifurcation of the elUptic fixed point, and 
it initially moves towards smaller k values. While it is traveling in the "wrong" direction, 
the elliptic 3/10 orbit has a winding number that does not exceed 1/5000; therefore orbits 
that bifurcate from it are of period 50000 or higher. We suspect that these would be 
orbits that are not smoothly connected to the AI limit, and therefore violate our working 
hypothesis. Since they bifurcate off the 3/10 orbit which in turn is smoothly connected, 
they are, however, at least continuously connected to the AI limit, so that the no-bubble 
conjecture could still be true. 
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In the dissipative case, the lowest period example we were able to find of 
a periodic orbit that is not smoothly connected occurs for b = —0.46, where 

at /c ~ 1.0346 the period 8 orbit with sequence (— ^ H has a period 

doubling bifurcation that creates a period 16 orbit that is not smoothly 
connected to the AI limit. Since this is the start of a period doubling 
sequence resulting in a what appears to be a strange attractor at a slightly 
smaller value of A;, we expect there are many saddle- node bifurcations in this 
region that create orbits as k decreases that are (presumably) not connected 
to the AI limit. This is reason to believe that the "no bubbles" conjecture 
is false when |6| < 1. 

In any case, continuation from the anti-integrablc limit has the advantage 
that at the beginning point all orbits exist, and they all continue nondegen- 
erately. 

4 Symbolic Dynamics 

In this section we introduce some notation for symbol sequences and bifur- 
cations. For simplicity we concentrate mostly on the area preserving case, 
b = I, though many results apply generally. 

Orbits in the anti-intcgrable limit are bi- infinite sequences s S S. When 
it is needed, we will indicate the current time along an orbit using a 

so that s = . . . S-2S-1.S0S1S2 The dynamics on s G S are given by the 

shift map, a : H ^ T, defined as 

Ct(. ..S_i.SoSlS2...) = . . . S_iS0"SlS2 • • • 

An orbit of the symbolic dynamics is periodic if the sequence s is periodic. 
We will denote an orbit of least period n by the string of n symbols and a 
superscript 00 to represent repetition: 

(sqSi . . . S„_i)°° = . . . Sn-2Sn-l ' SqSi ■ ■ ■ Sn-lSQ ■ ■ ■ 

Of course any cyclic permutation of a periodic orbit gives another point on 
the same orbit. 

Trivially, the map a has two fixed points, (+)°^ and (— )°°, and these 
correspond to the two fixed points of the Henon map. These are born in a 
saddle-node bifurcation at A; = —(1 -t- 6)^/4, which we denote by 

sn{(+r, (-)-}. 

We denote bifurcations with the general template 

parent — type {children} , 

where parent refers to the orbit that is undergoing the bifurcation, if any, 

and type is one of sn, pf, pd, or m/n, corresponding to a saddlc-nodc, 
pitchfork, period doubling, or rotational bifurcation, respectively. The set 
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of orbits created in the bifurcation is listed as the children. When the 
stability of these differ, we adopt the convention that the unstable child is 
listed first, and the stable one second. 

When b = 1 the fixed points of the Henon map are located at 



±Vl + ^ + e = e(l±VTTk 



z± 

The stability of a period n orbit of an area preserving map / is conveniently 
classified by the "residue" defined as 

R = \i2- TriDD) , 

so that an orbit is hyperbolic if it has negative residue, elliptic when < 
i? < 1 and is reflection hyperbolic when R > 1 |32|. The residues of the 
fixed points are 

R± = T^^l + e2 = T^Vl + k , (8) 

so the sign of the symbol is opposite to the sign of the residue of the fixed 
point. Thus the orbit (+)°° is always hyperbolic, while the orbit (— )°° is 
reflection hyperbolic for small e, or large k, but becomes elliptic at e = l/\/3, 
or A; = 3. 

For 5 = 1 the sequence (H — )°° corresponds to the period two orbit 



(+-r : (2o,^i) = (v/l-3e2-e,-Vl-3e2-e). 

This orbit exists only for e < 1/ \/3, and is created by a period doubling of 
the elliptic fixed point (when i?_ = 1). We denote this bifurcation by 

(-r-pd{(+-r}. 

Similarly there are two period three orbits. 



(- + +r : (zo,^l,^2) = (-^/^^-e,\/^^,^/^^). 

These are created in a saddle- node bifurcation at A; = e = 1; 

sn{(--+r, (- + +)-}. 

We list the low period orbits and their bifurcation values in Table ||. 
Another class of bifurcations shown in the table are rotational bifurcations. 
A rotational bifurcation occurs when the winding number of an elliptic orbit 
becomes co = m/n; we denote such bifurcations by the winding number of 
the parent orbit. For example the birth of orbits with winding number 1/n 
at the fixed point (— )°° is denoted 

(_)°o ^ 1/n {(- - +"-2)oo^ (_ + ^n-2)oo| _ (g) 



4 SYMBOLIC DYNAMICS 



10 



Parent 


Type 


Child 


Child 


A;- Value 




sn 


(_)oo 




-1 


i-r 


pd 


(+_)oo 


3 




sn 


(_ + _)oo 


(+ - +)~ 


1 


i-r 


1/3 


(_ + _)oo 




5 
4 


i-r 


1/4 


(+__+)oo 


(+ - ++r 





(+_)oo 


pd 




(- + — r 


4 


i-r 


1/5 


(- + + + -)°^ 


(+ + - + +)°° 


7-5^5 
8 


i-r 


2/5 


(-- + --r 


(- + - + -)°° 


7+5v^ 
8 




sn 


(+- + -+)°° 


(+ — +)°° 


5.5517014t 


i-r 


1/6 


(_+4_)oo 


(+ + -+^)°^ 


3 
4 


(+-+)°° 


pd 




(+_4+)oo 


5 
4 


(+ +)°° 


pf 


(+ + - + — )°° 


(- - + -++)°° 


3 




sn 


(-- + -+-)°^ 


(_ _+_3)oo 


3.7016569^ 


(+_)oo 


1/3 


(-- + -+-)°° 




15 
4 




sn 


(+_+3_)oo 


(_ _+3_)oo 


5.6793695^ 



1" 16P - 108A;^ + 105/c^ + 27A:^ - 97/c - 47 

1 16fc6 - 136fe5 + 213A:^ + 220P + 126A;2 + W8k + 81 



Table 1: Periodic orbits of the Henon map up to period 6 and their bifur- 
cations when 5 = 1. In the "type" column, "sn" indicates a saddle- node 
bifurcation, "pf a pitchfork bifurcation, and "pd" a period doubling bifur- 
cation. A rotational bifurcation is denoted by m/n, referring to the winding 
number of the parent at the bifurcation. For 1/3 the child is not created in 
the bifurcation, it exists before and after the bifurcation. If there are two 
children, the one listed in the first column has negative residue just after 
birth (except for the pf case). The real roots of the polynomials in the last 
rows give exact bifurcation values for the three approximations shown. 

This particular rotational bifurcation occurs when the multipliers of the 
fixed point are e*'^'^'^ or using Eq. (^), when k is given by 

= cos(27rw)(cos(27ru;) - 2) . (10) 

We have empirically identified the symbol sequences for rotational bifur- 
cations, and will present the general symbolic formula for these and for 
rotational "island around island" orbits in p9| . 

The residue of any periodic orbit of a Lagrangian system is easily com- 
puted from the matrix M formed from the second variation of the action 
[p3| . For a period n orbit of the Henon map this formula gives: 

where M is the periodic tridiagonal matrix with elements 



Mt,t-i = -be , Mt,t = 2zt{e) , Mt,t+i = -e 
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As we approach the anti-integrable hmit, z(e) — > s as e ^ and AI ap- 
proaches the diagonal matrix Diag(2sj). Thus we see that the residue be- 
comes infinite at the anti-integrable limit and its sign is given by — rir=o^ ^t- 
Hence, 

sign(i^(s)) = -(-l)^ (11) 
where j is the number of minus signs in the symbol sequence s. 

5 A Subshift of Finite Type 

In this section, we extend Theorem |^ to the case of a subshift of finite type. 
In particular, the biggest restriction in the proof of the theorem arises from 
the fact that the lower bound on the operator T given in Eq. (^) is weakest 
when the signs Sj+i = Si-i = —1 for positive b. We can improve the bound 
by restricting the set of admissible symbol sequences to forbid this particular 
case. The shift map restricted to this subspace is a subshift of finite type 
with the forbidden set = { — I — ,— — — }; that is, we define the shift 
space 

SjF = S\{s:3tGZ such that Sf_i = St+i = — } . 

This subshift can be easily described as a subshift on 2-blocks represented 

by { , — h, H — , ++}^. The subshift on the two-block space is represented 

by a vertex graph with the state transition matrix [^4| 



S = 



1 





1 




1 





1 














1 


V 


1 





1 / 



which indicates by Sij = 1 an allowed transitions from state j to state 
i. The two zeros and 8^2 come from the forbidden sequence in Sjf, the 
remaining are obtained because successive two-blocks overlap in one symbol, 
i.e., {st-ist) has to be followed by {stSt+i). The number of fixed points of 
period n for the subshift is given by 

Tr(5") = 7" + (1 - 7)" + 2(-l)"/2(j^ _ 1 mod 2) , 

where 7 = (1 + \/5)/2 is the golden mean. Thus the topological entropy for 
SjF is In 7. The number of distinct periodic orbits can be obtained from the 
trace formula by subtracting the number of periodic orbits for all factors of n 
and then dividing by the number of cyclic permutations, n. For comparison 
with the full shift and with the numerical results below, we give a list of 
these in Table |^. For example there are a total of 1, 465, 020 periodic points 
of the full shift with period n < 24, while there are only 12, 216 in the 
subshift SjF. 

When b is non-negative, orbits with symbol sequences in the subspace 
Tijr can be shown to persist longer than a general orbit: 
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l-'fjT'i r^/i 
1 (31 lUU. 






l-'oT'i r\ c\ 


y 




1 
1 


9 


1 




oou 




9 


1 
1 


n 

u 




1 1 fi1 


(JO 


3 


2 


1 


15 


2182 


90 


4 


3 


2 


16 


4080 


135 


5 


6 


2 


17 


7710 


210 


6 


9 


2 


18 


14532 


316 


7 


18 


4 


19 


27594 


492 


8 


30 


5 


20 


52377 


750 


9 


56 


8 


21 


99858 


1164 


10 


99 


11 


22 


190557 


1791 


11 


186 


18 


23 


364722 


2786 


12 


335 


25 


24 


698870 


4305 



Table 2: Number of orbits with minimal period n of the 2-shift and the 
subshift Ejr- 



Theorem 3 (Existence and Uniqueness of Ejr orbits). Suppose < 
b < 1. For every symbol sequence s G Sjf there exists a unique orbit z(e) of 
the Henon map Eq. ^) such that z(0) = s providing < e < emax; where 



-62 + 25 + 5 - 2V5T46 

[12) 



+ (i_5)(5_5) 

This theorem follows from the same argument that gave Theorem |l] with 
only minor modifications. We summarize the changes in the argument in 
the following discussion. 

Proof: When < 6 < 1, s E Sjp, and z G Bm, we can bound the norm 
of iterates of T in Eq. (^) using the inequalities 

ak < ||T'=(^)||oo </3fc , 
where the coefficients and Pk are determined by the recursions 

Pk = y/ l + e(l + b)P~i , 
Oik = \/l + e{ban-i - Pn-i) , 

with (3q = \ + M and uq = 1 — M. The sequence (3^ is the same as that in 
[11 1; it has the unique attracting fixed point 

/3oo = ^{e(l + 6) + vWW+4} . 

Since the recursion for {a„} depends on /3, but not vice versa, the coupled 
system also has a unique attracting fixed point, which is given by (aoo,/?oo) 
with 



i 1^6+^^252 + 4(1 _,^^)| 
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This imphes that for large enough n, T" maps the ball -Bi-q^o into itself. 
The map T is a contraction map on this ball providing ||DT'^||oo < 1- This 
leads to the same bound as that in |jll|, namely: 

(1 + 6) < 2aoo . 

After some silnplification, this inccjuality yields the formula, for c^iax* 

The 

condition that the operator T be real, e < l/(/3oo — baoo), is satisfied when- 
ever the map is a contraction. □ 

Similar arguments for negative b lead to bounds for a subshift with the 
forbidden subsequence + *— (where * is any symbol). However, this subshift 
is not of much interest, since there are only three periodic orbits in it. 



6 Horseshoe Boundary 

Theorem || provides an analytical bound on the parameter range for which 
the Henon map has a hyperbolic horseshoe. This bound corresponds to the 
dark grey region in Fig. According to Theorem |3| the subshift Sjf exists, 
in addition, in the lighter shaded area in the figure. This bound is valid only 
for 6 > 0, and meets the former at 6 = 0. 

Here, we use our continuation method to estimate the boundary of ex- 
istence of the horseshoe by following all orbits up to period 24 from the 
anti-integrable limit. In order for the numerical boundary to be valid we 
only need to assume that we can extrapolate from 24 to oo. Since there 



are at most 2" periodic points of period n in the Henon map |10|, we know 
that we are not missing any orbits because we have them all at the AI limit. 
This is no longer true after the first bifurcation because orbits that have 
disappeared might reappear for smaller k. 

To construct a numerical approximation for the boundaries, we first 
generate all symbol sequences for orbits of periods up to 24. Then, for fixed 
6, we numerically continue each orbit in e away from the anti-integrable limit 
and monitor its multipliers to detect bifurcations. For each b we record the 
smallest value of e at which a bifurcation occurs. The resulting numerical 
bounds in Fig. |l| are shown as solid curves. 

The numerical bound for the full shift is similar in shape to the ana- 
lytical one, but shifted to the right in e. While the analytical bound is 
symmetric under 6^—6, the numerical results are not. For example the 
first bifurcation at 6 = 1 occurs for e ~ 0.41888, while at 6 = —1 it occurs 
for e ^ 0.40167. In the logistic limit (6 = 0), Eq. (|l]) reduces to the logistic 
map, 

zt+i = -{z^-l), (13) 



for which the first bifurcation occurs at e = 1/^/2, where the orbit of the 
critical point becomes bounded. 
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Figure 1: First bifurcations for the Henon Map. The dark shaded region represents 
Theorem |l| and the hghter that of Theorem ^. The curves represent the numerical 
results for the first orbits destroyed up to period 24. Bounds for the subshift Sjr 
are indicated with a triangle symbol. 



When b is positive, the symbol sequences for the first pair of orbits 
destroyed up to period 24 extrapolate to orbits that are homoclinic to the 
fixed point (+)°°; we conjecture that these are the first orbits destroyed as 
e increases from 0: 

Conjecture 2. For positive b, the first bifurcation as e increases from 
corresponds to the homoclinic saddle-node bifurcation 

sn{+°°-(+)-+°°,+°°-(-)-+°°} . (14) 

The parenthesis in the middle enclose the "core" of the homoclinic or- 
bit, see the next section. A theorem of Smillie ||3^ implies that the first 
bifurcation destroying the Henon horseshoe must be a quadratic homoclinic 
tangency for some orbit. Our observations imply that it is a homoclinic 
bifurcation of (+)°°. When 6 < 0, however, the most natural description of 
the first bifurcation is as a heteroclinic tangency, which leads to 

Conjecture 3. For negative b, the first bifurcation as e increases from 
corresponds to the heteroclinic saddle-node bifurcation 

sn{-°° + (-) - +°°, + (+) - +^} . 

This does not contradict Smillie's theorem, as there are many homoclinic 
bifurcations that accumulate on this heteroclinic bifurcation. For example, 
for each m the orbits (— ™ + H |_™)oo — |_m-joo^ ^^^^ homoclinic 
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to the periodic orbit (— ™ + H |_m-joo^ ^^^^ bifurcation points of these 

homocHnic orbits hmit on that of the heteroclinic orbits as m — > cxd. 

Filtering the symbol sequences to choose only those in the subshift Sjr, 
we can use the same numerical data previously described to find the first 
bifurcation amongst the orbits in Sjr. This gives the solid curve marked 
with triangles in Fig. ||. This curve has a qualitatively different shape than 
the analytical bound. 

For reference we indicate in Fig. |l| the point k = 1.4, and b = —0.3, 
corresponding to the much studied Henon attractor. We also sketch the 
parameter range (b small enough, 1 < A; < 2) for which the theorem of 
Benedicks and Carleson implies that the Henon map has a transitive 
attractor with positive Lyapunov exponent. 

Note that the numerical horseshoe boundary does not depend on the no 
bubble conjecture. This is so because it is known that there can be at most 
2" periodic points of period n |l^]. Since we follow all of them up to period 
24 there can be no other orbits up to that period. In other words, orbits 
first have to be destroyed before they can be reborn. 

7 Homoclinic Orbits 

In this section we use the symbolic dynamics to classify orbits of the Henon 
map that are homoclinic to the hyperbolic fixed point p = (+)°° and study 
their bifurcations. We begin with some general terminology, referring to the 
Henon map as an example. 

Let / be an orientation preserving map^ of the plane with hyperbolic 
fixed point p. The stable and unstable manifolds of p are denoted by 
and W^, and a closed segment of such a manifold between two points a 
and P by Ty"[a,/3]. We use a parenthesis to denote an open endpoint of a 
segment. A segment that extends to the fixed point, e.g. W^{p,a], is called 
an initial segment of the manifold. The set of homoclinic orbits is the set of 
intersections D W^. A point a is on a primary (or principal) homoclinic 
orbit if the two initial segments to a touch only at a, i.e., 

W''{p,a]nW'{p,a] ={a} . 

Thus the initial segments to a primary homoclinic orbit define a Jordan 
curve; we call the interior of this curve a resonance zone. More generally, a 
resonance zone is a region bounded by alternating initial segments of stable 
and unstable manifolds [ 

For example, in Fig. |^ we sketch the left-going branches of the manifolds 
from p = (+)°° for the area preserving Henon map. There are precisely two 
primary homoclinic orbits; in the figure, we label points on these orbits with 
a and C- We choose to use a to construct the resonance zone. 

The orientation reversing case could be included by considering since its manifolds 
have the same geometry as those of /. 



37, 
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Figure 2: Stable and unstable manifolds for the Henon map at fc = 5 and b = 1 
shown in (z, z') coordinates. 

The stable manifold is divided into two invariant branches by the fixed 
point. An ordering is defined on each branch of W^, so that [3 <s "f for 
two points on a branch of if (3 is nearer to p along than 7, i.e., 
(3 G 7). We similarly define an ordering on each branch of W^. 

A segment of a manifold from a point to its iterate, W^{I3, /(/?)], is called 
a fundamental segment The union of the iterates of a fundamental 

segment is the entire branch of the manifold that contains (3. Moreover, 
since the iterates are disjoint, every homoclinic orbit on this branch must 
have precisely one point on each fundamental segment. 

For the Henon map, we focus on the left-going branches of and 
and the fundamental segments between f~^{C) and C- These also form the 
boundaries of the incoming and exit sets for the resonance zone defined by 
a [|3^. The exterior halves of these segments, W'^{aX) and H^"(/~-^(C), a), 
contain no homoclinic points since orbits on these segments are unbounded, 
so it is sufficient to look for homoclinic points on the interior halves, 

S ^ WV-\C),a]- (15) 



Every homoclinic orbit must have exactly one point on both S and U. 
Homoclinic orbits can be classified in a number of ways. The type |38], 
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of a homo clinic point /? i^ 

type(/?) = sup{j > : W'{p, n W^{p, /3] / 0} ; 

i.e., the number of iterates for which the stable initial segment to f^{(3) 
intersects with the unstable initial segment to (5. The type of a homoclinic 
point is invariant along its orbit. Primary homoclinic points have type 0. 

Homoclinic orbits on particular branches of and can also be clas- 
sified by their transition time. In general this is defined relative to a choice of 
a primary homoclinic point, C, and the fundamental segments W^"(/~^(C)) C] 
and C]- Any homoclinic orbit on these branches has exactly one 

point, /3, on the unstable segment. The transition time is the number of 
iterates required for /5 G 1^"(/~"'^(C)) C] to reach the stable segment: 

ttran.(/3) = A; if f\fi)eW\r\Q,Q] 

The value of the transition time depends upon the choice of fundamental 
segments, so it is not as basic a property as the type. 

In the simplest case, the transition time is easily related to the type of 
the orbit Q: 

Lemma 4. Assume there are exactly two primary homoclinic orbits, Q and 
a, and the segments S and U defined in Eq. ( fl^ j contain all of the homoclinic 
orbits. Then for each homoclinic point in (3 £ U, ttrans(/5) = type(/3). 

Proof: If /3 G [/ is of type i , then by definition W (p, /* {(3)] D {p, /3] ^ 
0. Now since a C and W''{p,C) D C) = 0> this implies that 

a <s However, n W{p,l3] = 0, which means that 

<s a, but there are no homoclinic points on W^{C,,a), so actually 
<s C- Now S contains every homoclinic point that reaches W^{p, C) 
in one iteration, so G S". □ 

Each homoclinic orbit has a Poincare signature that determines the di- 
rection of crossing of and at points on the orbit. We define the 
signature to be -|-1 if, looking along the unstable manifold in the direction 
of motion, the stable manifold crosses the unstable from the left to the right 
side. Crossings in the opposite direction have signature —1. If the manifolds 
do not cross but only touch (a topologically even intersection), the signature 
is defined to be zero. Since the map is orientation preserving, the signature 
is invariant along an orbit. Thus in Fig. |^ a and (" have signatures —1 and 
-|-1, respectively. The signature of a particular homoclinic orbit is typically 
not preserved in a bifurcation, but the total signature of the bifurcating 
orbits must be the same on each side of the bifurcation value. For example 

* Our definition of the type differs from Easton's slightly, to comply with his definition 
that type 1 is equivalent to the horseshoe. Rom-Kedar Ml uses the term Birkhoff signature 
instead of type. 
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a saddle-node bifurcation creates a zero signature orbit that splits into one 
positive and one negative signature orbit. 

For the Henon map, the AI symbol sequence can be used for the classifica- 
tion of homoclinic orbits. It is easy to construct homoclinic and heteroclinic 
orbits using the symbolic dynamics: an orbit heteroclinic from a periodic 
orbit (s)°° to a periodic orbit (s')°° has a symbol sequence that begins with 
a head sequence and ends with a tail sequence (s')°° with some arbi- 
trary, finite symbol sequence separating the head and tail. For example, the 
simplest orbits homoclinic to p = (+)°° are the primary homoclinic orbits: 

c = +~-.+~, 

a = +~ - . - , (16) 

corresponding to those we labeled in Fig. ^. These symbol sequences arise 
because as e — > the point a moves to the point — .— , while C moves to 
-.+ and /"HO to + .-. 

All other orbits homoclinic to p can be written in the form +°° — (s) — 
+°°, where s, the core, is any finite sequence — thus there is a one-to-one 
correspondence between finite symbol sequences and potential homoclinic 
orbits (all of which exist in the AI limit). This implies, for example, that 
near the anti-integrable limit there are 2^ homoclinic orbits with core length 
k. We will often denote an orbit homoclinic to p simply by writing the core 
in parenthesis, (s). Note that a given core (s) is not equivalent to any core 
with the same s cyclically permuted. 

The classification of homoclinic orbits by their symbol sequence can be 
used to compute other invariants. To determine the type of an orbit, we 
simply note that the AI symbols give exactly the same coding for an orbit 
as the standard symbolic coding for the horseshoe. This implies that the 
point +°° — . (s) — corresponds to a phase point on U, and the point 
+°° - (s) . - +°° is on S, thus 

Lemma 5. The transition time of the homoclinic orbit close to the AI limit 
is given by the length of the core sequence. 

For example, the homoclinic orbit +°° — ( h) — has the core sequence 

(— — h), and therefore has transition time 3. 

Similarly the signature of a homoclinic orbit in the horseshoe is given by 
simply counting the number of — signs in the core sequence. 

Lemma 6. The signature of a homoclinic orbit with core (s) close to the 
AI limit is given by —(—I)-' where j is the number of — signs in s. 

Thus the orbit (— — h) has signature —1.0 We will see that, when 6 = 1, 
some homoclinic orbits undergo pitchfork bifurcations, which change their 
signature, so this rule is not valid for all parameter values. 

^Eq. ( p^ ) implies that the signature is the same as the limiting sign of the residue of 
periodic orbits that approximate the homoclinic orbit. 
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Figure 3: Ordering of homoclinic orbits of type 1, 2 and 3 at fc = 6.25. The 
enlargement on the left shows the core sequences for the 14 homoclinic orbits up 
to type 3 The type orbits a and ( are not listed. 

The positions of the homoclinic orbits on U for orbits of type 1, 2 and 
3, labeled by their core sequences, are shown in Fig. |^ The order of the 
homoclinic orbits on the segments S and U close to the AI limit must be 
the same as the corresponding ordering of homoclinic points in the complete 
horseshoe by continuity. This ordering is equivalent to that of the logistic 
map, Eq. (|l3|), for the orbits forward asymptotic to the + fixed point. This 
gives an easy way to compute the ordering, see Fig. ^ In the logistic limit, 
all of the sequences forward asymptotic to the fixed point are destroyed 
when the orbit of the critical point becomes bounded at e = 1/V2. The 
fixed point has a single preimage, which is a = . — Every other 

orbit that is forward asymptotic to .+°° has the form .(s) — \-°°. 

In the area preserving case the ordering of the orbits along S is equivalent 
to that on U upon time reversal. Thus a type t point +°° — (siS2 . . . st) . — \-°° 
on S is in the same relative position as the point +°°— . (stSt-i ■ ■ ■ si) — on 
U. Close to the AI limit we always have this ordering on the manifolds, which 
is just another way of saying that the map is conjugate to the horseshoe map. 

So long as there are no homoclinic bifurcations, then the orderings 
and >s are just given by the usual unimodal ordering as stated in 

Lemma 7. The ordering on U and >s on S close to the AI limit is 
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a 



-++ 



.+-+ 



.++- 



.+++ 



Figure 4: Ordering of pre-periodic points for the + fixed point of the logistic limit 
of the Henon map. The symbols are determined by the itinerary of the orbit relative 
to the critical point at z = 0. 

given by 



.e + ... 



.o + ... 



+ e. 

+ 0. 



+ 



oo 
oo 



. e — 
. o — 



e . 
o . 



where e / o are finite sequences with an even / odd number of minus signs, 
respectively. 

The ordering shown in Fig. |3| and Fig. Q is exactly this one upon ap- 
pending the "homoclinic tail" — \-°° to the cores. The maximal orbit on U 
is corresponding to the tail .+°°, the minimal orbit is a, corresponding 
to . - +°°. 



8 Homoclinic Bifurcations 

Homoclinic bifurcations are bifurcations between homoclinic orbits. Com- 
pared to ordinary bifurcations of periodic orbits they possess additional 
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structure because the invariant manifolds (with their ordering) must be in- 
volved in the bifurcation process. To make this explicit, we say that two 
homoclinic orbits (3 and 7 are double neighbors if the segments TF"[/9,7] 
and Ty*[/3, 7] contain no other homoclinic orbits. Three ordered homoclinic 
points (3 <u ^ <u S are triple neighbors if both /?, 7 and 7, 5 are double neigh- 
bors. An obvious observation with nevertheless important consequences is 
the "double neighbor" lemma: 

Lemma 8. Two homoclinic orbits [3 and 7 cannot bifurcate unless they are 
double neighbors. 

The converse gives a simple forcing relation: before (5 and 7 can bi- 
furcate any homoclinic orbit on either segment between them must have 
disappeared. 

Another consequence is the transition time lemma: 

Lemma 9. If two homoclinic orbits (3 and 7 bifurcate then they must have 
the same transition time itrans- 

Proof: Let (5 and 7 be neighbors on [/. If their transition time is dif- 
ferent then they are not neighbors on S*, so they cannot bifurcate. □ 

This allows us to extend Lemma ^ away from the AI limit, so that one 
can take the transition time as an adequate replacement of the period: 

Corollary 10. The transition time of a homoclinic orbit never changes. 

Proof: Since the transition time is an integer it cannot change under 
smooth deformations. It could only change at bifurcations, but we have just 
seen that only orbits with the same transition time bifurcate. □ 

Therefore homoclinic bifurcations only take place between double neigh- 
bors with the same transition time, i.e., core length. Close to the AI limit 
the horseshoe is still complete. In this situation it is possible to find all 
neighbors: 

Lemma 11. Two homoclinic orbits on U are neighbors in the complete 
horseshoe if and only if they are of the form +°° — . (s+) — and +°° — 
.{s-)-+°°. 

Proof: We have to show that there is no homoclinic orbit with core 5 
such that (0+) <u {S) <« (o— ) or (e— ) <„ (5) (e+), where e = s if s has 
an even number of minus signs or o = s if this number is odd. If the initial 
string in 6 differs from s then 6 can not be between the sequences (s— ) and 
{s+), therefore 6 = s . . . . It is simple to see that 

• 6 ~t~ * * * ' e -\- h , • 6 • • • • ^ f~ ? 
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and similarly 

.oH <u-o + -+°°, .0 >«.o--+°°. 

Since 5 = s . . . it must be of one of the forms on the left hand sides, but 
then the inequalities show that it is not between (s+) and (s— ) hence they 
must be neighbors. Conversely, suppose we have two neighboring homoclinic 
orbits (on U) .a and .b with .a <u They must differ in at least one 
symbol so call the first such difference x. Their leading common symbols 
are denoted by s, so that a = sxa and b = sxj3 for some sequences a and /?, 
where x is the opposite symbol to x. Applying the ordering relation to the 
possible combinations of s and x gives either 

.e- a <u . {ey) - +°° <u •e + [3 or . o + a <„ . {oy) - +°° <„ .o- P , 

where the choice of the symbol y depends on whether s is even or odd and 
whether a and /3 are — • Specifically, choose y = + if either s = e and 
/? 7^ — \-°° or s = o and a / — \-°°. Choose y = — if either s = e and 
a 7^ — -)-°° or s = o and /? 7^ — +°°. If neither a nor /? are — +°° either 
choice for y works. When either a, /?, or both differ from — we have 
constructed an orbit . (sy) — which is between . a and . b — hence . a and 
.b are not neighbors. But this is a contradiction so a, P = — \-°°. □ 

For a bifurcation to occur it is not enough that the orbits be neighbors 
on U, but they must be double neighbors. In the reversible case this almost 
gives the proof of Conjecture ^, but here we are working in the smaller class 
of orbits homoclinic to p, the hyperbolic fixed point. So far we did not make 
use of the reversibility of the map, i.e., the results are valid for all b. From 
this point on we will always only talk about the area preserving case. Note 
that the ordering is used in a range of parameters before the first bifurcation 
occurs, so the horseshoe ordering is still valid. 

Theorem 12. In the area preserving Henon map the first homoclinic bifur- 
cation of the invariant manifolds of the fixed point (+)°° is 

sn{+°°-(+)-+~,+°°-(-)-+°°} . 

Proof: By Lemma |ll| we know that all neighbors on U in the complete 
horseshoe are of the form (sit). For these sequences to be double neighbors 
they must be neighbors on S as well. By reversibility this is equivalent to 
the sequence and its reverse being neighbors on U. But this only true if s is 
empty. The only double neighbors in the complete horseshoe are therefore 
the two orbits +°° - (+) - +°° and +°° - (-) - +°°. Therefore they must 
bifurcate first. □ 

To approximate a homoclinic orbit, which possesses an infinite number of 
points in phase space by a periodic orbit with only a finite number of points 
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Table 3: Bifurcations in periodic approximations to the homoclinic type 1 
orbit, which is the first orbit destroyed for 6 = 1. Here we use a * to denote 
both + and — , giving both orbits involved in the bifurcation. 



we require that the Hausdorff distance of these two point sets vanishes as the 
period approaches infinity. Thus for an orbit homoclinic to (+)°°, we study 
a sequence of approximating periodic orbits with an increasingly long string 
of + symbols. In particular the rotational orbits given in Eq. (^) converge 
to C and a in the limit. 

In Table ^ we list the first 11 members of the sequence approximating the 
transit time 1 homoclinic orbit, and the corresponding sequence of values, esn 
at which these orbits undergo a saddle- node bifurcation when 6 = 1. These 
values converge geometrically to the parameter at which the homoclinic or- 
bits bifurcate, and the ratio of successive differences ( a "Feigenbaum ratio" ) 
is computed in the fourth column of the table. As is known theoretically for 
6 < 1 pT| , the convergence rate, 5, approaches A, the multiplier of the 
fixed point p. From our data, the convergence rate 6 agrees up to 6 digits 
with the multiplier 

A K 7.0344478 



of the fixed point p when k ~ 5.699310786700. Thus, our observations 
indicate that the convergence rate is given by the multiplier in the area 
preserving case as well, where to our knowledge no proof exists. 

The third column in the table is the extrapolation for the converged k 
value, given by Aitken's method 

k -k 



A^ik„ 



where A is the forward discrete difference operator. Thus we see that there 
is a saddle-node bifurcation of the type 1 homoclinic orbits, 

sn{+°°-(+)-+°°,+°°-(-)-+°°} , 
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at 

esn(l) ~ 0.418879233367 or A;sn(l) ~ 5.699310786700 . 

This also corresponds to the parameter value at which the topological horse- 
shoe for the Henon map is destroyed, and is the value in Fig. [l| at b = 1. 

Since there is a sequence of saddle-node bifurcations that limit on the 
homoclinic bifurcation, there are elliptic islands arbitrarily close to the de- 
struction of the horseshoe. This corresponds to an area preserving version 
of the results of Gavrilov and Silnikov |20, ^] and Newhouse [l^ . 



Since we can in principle follow every finite orbit from the anti-integrable 
limit we can begin to study the sequence of bifurcations that occur after the 
horseshoe is destroyed, see Table ^. For example, the bifurcation diagram 
for all of the homoclinic orbits of type three or less is sketched in Fig. ^. 
The vertical ordering in this sketch is the same as that on the segment U 
with a and C shown. The bifurcation diagram is highly influenced by the 
time-reversal symmetry of the area preserving Henon map — we will discuss 



this symmetry in the next section. As expected from the general theory |43], 
we observe three kinds of bifurcations: 



Symmetric saddle-node bifurcations resulting in the creation of a pair 
of type t homoclinic orbits with opposite signatures. For example, in 
Fig. |5|, the type 3 orbits with cores (+ + +) and ( — I — ) are born in 
such a saddle-node at /c 0.386. 

Pitchfork bifurcations of type t symmetric homoclinic orbits, creating a 
pair of type t asymmetric orbits that are related by time reversal. For 
example, the ( — I — ) orbit pitchforks at ~ 0.720 creating the orbits 
( — and (+H — ). A pitchfork bifurcation requires triple neighbors 
to occur. The parent orbit of a homoclinic pitchfork bifurcations is 
always created in a symmetric saddle-node bifurcation. Up to type 
11 there are are only 9 symmetric saddle- node bifurcations which do 
not undergo a homoclinic pitchfork bifurcation on their way to the AI 
limit. 

Asymmetric saddle-node bifurcations creating two symmetry related pairs 
of asymmetric orbits. This bifurcation first occurs at type 4. For ex- 
ample the two pairs {( — I ), ( — I h)} and {( 1 — ), (H 1 — )} 

are created at A; ~ 5.18. Generically, asymmetric saddle-node bifur- 
cations require two pairs of double neighbors to occur because of the 
symmetry. 

The shaded region in Fig. |5| represents the range of k for which the area 
preserving Henon map exhibits a horseshoe. Along the left edge we label 
each orbit with its core symbol sequence. 

The first type t homoclinic orbits are created by a saddle-node bifurca- 
tion when the segment f~*{S) first intersects U. We denote this parameter 
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Figure 5: Sketch of bifurcations in the homoclinic orbits up to type 3 (6 



Parent 


Type 


Child 


Child 


/s- Value 




sn 


(- + +-) 


(+ + ++) 


-0.133474 


(- + +-) 


pf 


(- + ++) 


(+ + +-) 


-0.044273 




sn 


(- + -) 


(+ + +) 


0.385556 


(- + -) 


pf 


(- + +) 


(+ + -) 


0.719630 




sn 


(— ) 


(++) 


1.627779 


(— ) 


pf 


(-+) 


(+-) 


3.091505 




sn 


(+ - -+) 


( ) 


3.98213640 


( ) 


pf 


( +) 


(+-— ) 


3.98213641 




sn 


(+-+) 


( ) 


4.706399 


(+-+) 


pf 


(--+) 


(+--) 


4.816792 




asn 


(- + -*) 


(* - +-) 


5.188561 




asn 


(*-++) 


(+ + -*) 


5.619922 




sn 


(+) 


(-) 


5.699311 



Table 4: Homoclinic bifurcations up to core length 4. 
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value by ksn{t)- This marks the creation of the subset of the incoming lobe 
of the turnstile with transition time t We observe that when 6 = 1, 

this homoclinic saddle-node bifurcation is 

sn{(+*),(-+*-2-)} atfesn(t). 

Following this, the orbit ( — — ) undergoes a homoclinic pitchfork bifur- 
cation at kpf{t), creating the pair 

(- -) - Pf {-+'-')} at fepf (t) . 

However, when b ^ 1, the initial symmetric bifurcation and the following 
symmetry breaking pitchfork are replaced by a pair of nonsymmetric saddle- 
node bifurcations. In this case the first type t bifurcation is the homoclinic 
saddle-node 

sn 

According to the double neighbor lemma, certain bifurcations cannot 
occur prior to other homoclinic bifurcations because the corresponding se- 
quences block other sequences from being neighbors. In order to determine 
which orbits are neighbors even beyond the first bifurcation we make the 
assumption that the following symbolic ordering conjecture holds: 

Conjecture 4. The symbolic horseshoe ordering on the invariant manifolds 
given in Lemma persists. 

The ordering relations give a unique construction of the order of the 
points on U and S, and this implies that a schematic construction of the 
intersections of /"*(S') with U can be constructed solely from a list of which 
orbits exist at a given parameter value. Such a schematic manifold plot is 
shown in Fig. ^, for all homoclinic orbits that exist at k = 5.53 up to type 
5. 

We can also construct a schematic bifurcation diagram for homoclinic 
orbits, as in Fig. ^, by drawing a horizontal line from /c = cxd to the k- 
value at which a particular homoclinic orbit is destroyed — actually we stop 
the figure at /c = 6, since there are no bifurcations for larger fc-values. We 
order the homoclinic orbits vertically according to their unimodal ordering 
on U as usual. In this bifurcation diagram the vertical connections indicate 
which orbits eventually do become neighbors and bifurcate. So as to avoid 
artificially crossing lines, we connect pairs of asymmetric saddle-nodes by 
lines at the right edge of the figure to indicate that they must bifurcate at 
the same fc-value. 

We say that a bifurcation straddles the centerline if the pair of orbits 
involved are on either side of center of the U ordering, or if one of the two 
pairs of an asymmetric saddle-node straddles the center line. 

Through type 6, each symmetric saddle node is followed by a pitchfork 
bifurcation, just as we observed in Fig. ||, with the exception of the very first 
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bifurcation, sn{(+), (— )}, which corresponds to the smahest loop straddhng 
the center in the figure. That this is in fact the smallest loop and therefore 
the first bifurcation is the content of Theorem |l^. 

Moreover, it is remarkable, but perhaps misleading, that through type 6 
every bifurcation straddles the center. Therefore all homoclinic bifurcations 
up to type 6 are forced by nesting around the center. In particular this 
means that their unimodal ordering gives the bifurcation ordering, like in 
unimodal maps. 

This simple forcing relation is destroyed with the appearance of a sym- 
metric saddle-node without pitchfork of type 7 (see Table |6[ Also at type 7, 
there is an asymmetric saddle-node quadruple which does not straddle the 
center. Interestingly enough, this is the same bifurcation that marks the 
upper k endpoint of one of the gaps that we discuss in § |lO| . 

It is difficult to visualize the full homoclinic bifurcation diagram for 
larger type orbits. To do so, we plot only the horizontal lines, to indicate 
the range of existence of an orbit; this diagram up to transition time 11 is 
given in Fig. ^. The approximate self-similarity in this picture seems to be 
related to some of the gaps we discuss in § |To| , namely those that are related 
to symmetric saddle-nodes without accompanying pitchforks of type 7, 9 
and 11. 
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Figure 6: Schematic drawing of U (dashed hne) and / '(S*) (solid line) up to type 
5 for fc = 5.53. 
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Figure 7: Bifurcation diagram of homoclinic orbits up to type 5 {b — 1). Types 
1,2,3 are shown as dotted lines (recall Fig. ||); type 4 is dashed; and type 5 is solid. 
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Figure 8: Existence plot of homoclinic orbits up to type 11. For each homoclinic 
orbit a line is drawn from large k to the parameter value where this orbit is de- 
stroyed. The vertical position of each line is its formal position on U according to 
the unimodal ordering. 



9 SYMMETRIC HOMOCLINIC BIFURCATIONS 



30 



9 Symmetric homoclinic bifurcations 

As we mentioned above, the bifurcation diagram of the area preserving 
Henon map is restricted by the fact that the map has a time-reversal sym- 
metry. Here we briefly recall a few well known facts about reversible maps 
[44 1, and apply them to the study of homoclinic bifurcations. 

A map / has a time-reversal symmetry when it is diffeomorphic to its 
inverse by: 

Rf = f-'R . 

We call the map R a reversor for /. Often, as in our case, the reversor is 
an involution, R^ = I. Note that each of the maps f^R is also a reversor, in 
particular, we call fR the complementary reversor to R. The fixed set of a 
reversor 

fix(i?) = {x : Rx = x} , 

is of particular interest. For the case when R is an orientation reversing 
involution of the plane fix(i?) is always a curve that goes through infinity, 



thus dividing the plane into two pieces |45|. 

A reversor maps an orbit . . . Zt-i, Zt, Zt+i ... of the map onto another 
orbit . . . Rzt+i, Rzt, Rzt^i .... A symmetric orbit is defined as one that is 
mapped onto itself by R. It is easy to see that any symmetric orbit must 
have points on fix(i?) U fix(/i?) and conversely. Moreover, if the orbit is not 
periodic, it has a unique point on one of these fixed sets, and if it is periodic 



it has exactly two points on the fixed sets [46|. 

Reversible maps need not be area preserving, though the multipliers 
of an orbit and its symmetric partner must be reciprocals of one another. 
Application of this to the fixed points gives that the Henon map is reversible 
only when 6 = ±1. For a symmetric orbit reversibility implies that the 
product of the multipliers must be one. For the case 6 = ±1 a reversor is 
R{x,y) = {—y,—x), and a complementary reversor fR{x,y) = {—x — k + 
y'^iV). The fixed curves are 

fix(i?) = {{x,y) : X = -y] , 

MfR) = {(x,y):x = \{y^-k)}. 

Suppose that p is a symmetric, hyperbolic fixed point of a reversible map. 
Then, as pointed out by Devaney |^^, the stable and unstable manifolds of 
the map are related by R: 

Lemma 13. Let and he the stable and unstable manifolds of a 
symmetric fixed point p. Then RW'^{p, /3] = W^{p, Rj3\. 

Proof: By definition, when (3 G M^", then f~^{f3) p as t ^ oo. Then 
Rf-\(3) = f\R(i) Rp = p. Thus, Rf3 G W. Since R is a diffeomor- 
phism, RW''{p,/3] = W'{p,R(3]. □ 
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Corollary 14. // intersects the fixed set of a reversor, then the inter- 
section point is homoclinic. 

Homoclinic orbits of symmetric periodic orbits either come in symmetric 
pairs or are symmetric, and there must exist symmetric homochnic orbits: 

Lemma 15. Let p he a symmetric, hyperbolic fixed point, and [3 a homo- 
clinic point, and suppose that R is an orientation reversing involution. Then 
Rf3 is also a homoclinic point. Moreover, there exist symmetric homoclinic 
points on fix(i?) and fix(/i?). 



Proof: By Lemma y, since /3 G VK' n VK" then R(3 e W r\W' , so it 
is homochnic as weh. Since fix(i?) divides the plane and (5 and Rf3 are on 
opposite sides of this curve, the segment W'^[f5,Rl3\ must cross fix(i?), and 
the crossing point is necessarily homoclinic and symmetric. We can argue 
similarly for fR. □ 

As is well known, pitchfork bifurcations occur with codimension one in 
maps with a symmetry E^. This occurs for homoclinic bifurcations as well, 



as was suggested in [48|. We observed such pitchfork bifurcations in Fig. g. 
A pitchfork typically occurs after a symmetric, type f > 1, saddle- node 
bifurcation creates a "tip" of inside the entry lobe of the turnstile. As 
this tip grows, one would normally expect it to bend around, as sketched in 
Fig. ^, creating more type t homoclinic points by saddle-node bifurcation. 
In fact, it is a simple consequence of the linear ordering along and 
combined with reversibility that a single saddle-node bifurcation like that 
sketched in Fig. ^ is impossible: 

Theorem 16. Suppose that f is an orientation preserving, reversible map, 
with a symmetric fixed point p, and S and U = RS are segments of its 
stable and unstable manifold bounded by adjacent primary homoclinic orbits. 
Suppose that a pair of symmetric homoclinic points on U , (3 <s j are created 
in a saddle-node bifurcation. Then it is impossible for there to he a single 
saddle-node bifurcation as a tangency of W^{I3,^) with either piece of U \ 

Proof: Since /3 <s 7, and RW^ = W^, we have R(3 <„ Rj. Suppose 
that P and 7 have transition time t. Then € S, but since /? is symmet- 

ric this point must be the same as Thus f~^{R(3) = (3, and similarly for 
7. Since the ordering is preserved by iteration, then /3 <u 7. Now suppose 
there is a tangency at a point 6 = W^{f3, 7)n(C/\VF"[/3, 7]), i.e., not between 
P and 7. We sketch such a configuration in Fig. Thus (3 <s S <s 7. By 
symmetry, RP R6 Rj. Since the ordering is preserved by iteration, 
we have /3 <„ f~\R5) <u 7. Thus f-\R5) G VF"(/3,7) and so this point 
is not 5 (consequently the orbit of 6 is not symmetric). Since the manifolds 
are tangent at 5, they are also tangent at R5 by symmetry. Thus there is 
a second, simultaneous tangency, on U at f~^{R5) which contradicts the 
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fix(R) 




Figure 9: Impossibility of the second symmetric bifurcation as described in The- 
orem ^ Stable manifolds are shown as solid and unstable manifolds as dashed 
curves. Reflection of a tangency at S gives a tangency at RS that is ordered incor- 
rectly. 



assumption that a single tangency occurs. □ 



There are three possible resolutions: first one of the two orbits, /3 or 7 
could undergo a pitchfork bifurcation creating a symmetry related pair of 



homoclinic orbits. For example, Fig. |10| shows part of the homoclinic tangle 
at a parameter value where the type two homoclinic orbit with core sequence 

( ) pitchforks. As k increases this results in the creation of a pair of type 

2 orbits with cores ( — h) and (-1 — ), see Fig. |Tl|. Note that the new orbits 
are not symmetric, but that the reversal of (— +) is (+— ), so they form a 
symmetric pair. 

The second possible bifurcation is a single-saddle node on the segment 
W^{(3,^); this happens, for example, whenever a "tip" of an iterate of S 
returns to U. This first occurs at type 3, for the bifurcation sn{(* — *)}. 

We sketch a similar case, at type 4, sn{(* *)}, in Fig. 12 which occurs 

at 3.982. 

The third possible bifurcation is a pair of asymmetric saddle-node bi- 
furcations. This first occurs for homoclinic orbits of type 4. For example, 

the bifurcations sn{(-| 1 — ),( 1 — )} and its time-reverse, sn{( — h 

(- H )} occur at k ^ 5.1886. We sketch U and f~^{S) at this bi- 
furcation in Fig. [1^. This bifurcation also corresponds to the lower endpoint 
of an apparently hyperbolic parameter interval for the Henon map, as we 
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Figure 10: Stable and unstable manifolds for the (+)°° fixed point of the Hcnon 
map at & = 1 and k = 3.09151, where there is a cubic tangency of the manifolds at 
the +°° — ( ) — homoclinic orbit. 




Figure 11: Type 2 homoclinic orbits of the Henon Map at fc = 3.5. 
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■(-+3) 

Figure 12: Sketch of two possible homoclinic saddle-node bifurcations of type four. 

A symmetric saddle-node creating (* *) occurs on W'^{{ — h +— ), (+ + ++)) in 

(a). An asymmetric saddle-node creating (* — h— ) and ( — I — *) occurs with one 
point on (+ - -+)) in (b). 



discuss in the next section. 

Note that the antimonotonic bifurcations shown to exist in the area 
contracting case |13| are exactly forbidden by this theorem. 



A symmetric saddle-node followed by a pitchfork is a common bifurca- 
tion. For example, the parameter values, ksnit), at which the first type t 
homoclinic orbit is created decrease monotonically with t. Thus at ksn{t) 
the first type t orbit is born and there are no homoclinic orbits with type 
less than t. For t > 1, at ksnit — 1) the segment f^~^{U) must intersect with 
5, so that f^{U) intersects with f{S). In order for this to happen (when 
6 = 1), there is a pitchfork bifurcation for kpi[t) £ [ksn{t), kg^it — 1)] of the 
type t homoclinic orbit -|-°° — ( — — ) — giving rise to the pair of 
orbits with symbol sequences 

(-+*-2-)-pf {(-+*-!), (+*-^-)}. 

We see that the children of this bifurcation differ from their parent in a 
single symbol and they differ from each other in two symbols. Table |5| lists 
the first few such homoclinic bifurcation values obtained by extrapolation 
of the first few members of the approximating orbit sequence. 

The distance (in k) between the birth of the type t orbit and its pitchfork 
bifurcation shrinks to zero as the type increases. 
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I-^tIt' h Tr^iT" 1^ ( hilriT'tiTi 




1 


(*) 


5.69931078670 






2 


(**) 


1.62777931098 


(-+),(+-) 


3.09150542113 


3 


(* + *) 


0.38555621701 


(-+'),(+'-) 


0.71963023592 


4 


(* +^ *) 


-0.13347378530 




-0.04427324816 


5 


(* +^ *) 


-0.39678970175 


(-+'),(+'-) 


-0.36787481134 


6 


(* +^ *) 


-0.54918558488 


(-+'),(+'-) 


-0.53740149261 


7 


(* +^ *) 


-0.64623270965 


(-+'),(+'-) 


-0.64032496327 


8 


(* +6 *) 


-0.71262572399 


(-+'),(+'-) 


-0.70916824264 


9 


(* +^*) 


-0.76055766670 


(-+'),(+'-) 


-0.75830622014 


10 


(* +^ *) 


-0.79659407362 


(-+'),(+'-) 


-0.79501732767 



Table 5: Pitchfork bifurcations from the first type t orbits up to type 10. 



10 Intervals with no Bifurcations 

Davis, MacKay and Sannami (DMS) Q used the numerical method of Biham 
and Wenzel to compute the periodic orbits for the area preserving Henon 
map. They showed that up to period 20, there are intervals of parameter 
where there appear to be no orbits created or destroyed. They studied 
a particular parameter interval near the destruction of the horseshoe, and 
elucidated the symbolic dynamics of the corresponding homoclinic tangle. 
We will refer to this interval as the DMS gap. Though the method of Biham 
and Wenzel is guaranteed to work close enough to the AI limit [|lT[, it can 
fail [0. We tested the DMS results using our continuation technique. The 
use of parallel computation allowed us to extend the original experiment by 
an order of magnitude in size so that we followed all orbits up to period 
24 — recall from Table |2| that there are a total of 1,465,020 possible orbits. 
We verify the DMS results and identify the symbol sequences of the orbits 
that form the boundaries of the DMS gap. 

In our experiment we follow all orbits up to period 24 and record the 
minimal parameter values at which they are destroyed. We then assume 
that each orbit exists only up to that value of k; this procedure is not 
entirely correct, because a few orbits loop back and forth in parameter under 
continuation. This is related to the vanishing of twist in the neighborhood 
of a period tripling bifurcation [^]. However, the number of low period 
orbits for which this happens is very small. 



In Fig. 13 we show the number of orbits that exist as a function of k, with 
the caveat that no value is plotted if the number of orbits does not change 
from the previously plotted point. This plot is equivalent to that of DMS, 
except that we leave gaps in the intervals where there are no bifurcations. 

At the anti-integrable limit the map exhibits a horseshoe so all of the 
periodic orbits are present. As we move away from the anti-integrable limit 
we see a decline in the number of periodic orbits as orbits collide and are 



destroyed. Flat intervals in Fig. 13 represent intervals of parameter where 
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Figure 13: Number of periodic orbits of the Henon map up to period 24, for 
6=1. The maximal number is reached at A: « 5.69931078745, when the horseshoe 
is formed. The endpoints of the largest gap studied by DMS are labeled by L and 
R. 



very few bifurcations occur. Gaps in the plot indicate intervals of parameter 
where there are no bifurcations. The creation of the first type t homoclinic 
orbits gives rise to flat intervals. We observe that the left endpoint of each 
of the larger flat intervals for /c < 3 corresponds to /csn(i) for the saddle- node 
bifurcation of the first type t homoclinic orbits; these are marked in Fig. ^ 
and in the enlargement. Fig. 14. Similarly, the parameter values kpf(t) are 
also marked; note that these pitchfork bifurcations are located well beyond 
the right endpoints of the flat intervals. Each of the gaps in the flat intervals 
for k < 3 must eventually fill in if we go to high enough period because in 
this range of k the area preserving Henon map has an elliptic fixed point. 
Recall that an m/n bifurcation from the elliptic fixed point occurs at the 
parameter values k^/^ given in Eq. (|^), and these values are dense in the 
interval —l<k<3. Moreover, invariant circles bifurcate from the elliptic 
fixed point for each k^j for sufficiently irrational to. The same argument can 
be used up to the end of the period doubling cascade of the fixed point at 
k 4.13616680392, since each period doubling creates an elliptic orbit. 
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Left Endpoint Core 


kL 


Right Endpoint Core 


kR 


(- + -*- + -) 


4.55931896797 


(+ + -*-*- + +) 


4.59567964802 


{+''- + - *) 


4.84317164217 


(+^ -*- + - + + -) 


4.86795762007 


(- + -*) 


5.18851121215 


(+ + -*- + +) 


5.53765692812 


(- + + -*- + + -) 


5.56490867348 


(+ + -*-+^) 


5.60872105039 


(- + + -*) 


5.63190980280 


(+3 _ * _ +3) 


5.67769222229 



Table 6: Homoclinic bifurcations bounding the gaps in Fig. 13 



There are a number of distinct gaps in Fig. |l^; the 3 larger gaps were 
studied by DMS, especially the largest one, near k = 5.5 indicated by L 
and R in Fig. |l^. DMS conjecture that the dynamics in each gap is hyper- 
bolic, and consequently there are no bifurcations in a gap. Our numerical 
evidence, which extends their study by an order of magnitude, supports this 
conjecture. Upon examining the orbits that limit on the endpoints of the 
gap up to period 24, we can extrapolate and find that each of the five largest 
gaps is bounded by a homoclinic bifurcation, see Table |6[ Thus we see that 
the gaps do not fill-in with orbits converging on the homoclinic bifurcations, 
but we cannot rule out that there are other, unrelated period orbits with 
period larger than 24 that are created at parameter values in the middle of 
a gap. 

We observe that there are two types of bifurcations bounding the gaps: 
symmetric and asymmetric saddle-node bifurcations. The asymmetric saddle- 
nodes result in the creation of two pairs of homoclinic orbits, the one listed in 
the table, and its time-reverse. Typically we observe that symmetric saddle 
node bifurcations in homoclinic orbits are followed by pitchfork bifurcations. 
In fact we observe that among all of the homoclinic orbits through type 11, 
there are only 9 special saddle-node bifurcations which are not followed by 
a pitchfork bifurcation. We believe that each of these special bifurcations 
corresponds to the endpoint of a gap. For example, the first type 1 bifur- 
cation is not followed by a pitchfork, and it gives the left endpoint of the 
gap corresponding to the horseshoe. The four gap endpoints in Table ^ that 
correspond to symmetric saddle nodes are each of this special type. The 
four remaining special pairs are each of type 11, and of these at least two 
bound gaps of widths Ak ~ 6(10)~^. With our resolution it is not possible 
to clearly identify the final two as gap boundaries. 

In Fig. 14 we show an enlargement of Fig. 13 but also include the data 
from the subshift, Sjr. In the upper right corner of Fig. 14 we see the tail 
end of the exit time 2 plateau. We also labeled the first large gap in Sjf 
after the subshift is destroyed — this is the subshift analog of the DMS gap. 

As in the DMS gap, the left and right boundaries, denoted L and R, 
correspond to a pair of homoclinic saddle-node bifurcations with the core 
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Figure 14: Enlargement: number of periodic orbits up to period 24 for the full 
shift and the subshift Sjr. The left and right ends of each plateau correspond to 
fcs„(t), and kp{{t). 



sequences 

asn {(* - 
sn{(+' 



- + +-),(- + + -+*)} at A:sn(i^) ~ 1.533898312 , 
+ + at A:sn(-R) « 1.583387630 . 



Note that the right endpoint of the gap corresponds to an orbit whose part- 
ner is not in the subshift! 

The curves for all orbits and for the subshift are remarkably similar 
and it appears that the growth of orbits in the subshift gives an accurate 
representation of the overall growth of orbits in the full shift for this range 
of parameters. This is especially remarkable given that when all the orbits 
exist, the subshift contains less than 1% of the orbits in the full shift up to 
period 24. The figure shows that for small k, the number of orbits in the 
full shift is nearly a constant multiple of that in the subshift. 

Observing that the gaps are bounded by homoclinic orbits, we regener- 
ated the orbit growth plot using only homoclinic orbits. As expected, the 



gap structure and overall shape of Fig. |13| is almost completely captured by 
the homoclinic orbits alone. 
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11 Conclusions 

Continuation from an anti-integrable limit is an effective technique for study- 
ing orbits providing that there are no isolated bubbles in the bifurcation 
diagram. In ||Tl[| we applied the anti-integrable theory to the Henon map 
to obtain a new proof of the well-known analytical bound of Devaney and 
Nitecki In Theorem ^ we apply a similar argument to a restricted set 
of orbits to find an analytical bound for the existence of a subshift of finite 
type. We present both analytical bounds together with the optimal bounds 
generated numerically using our continuation method. We observe that the 
horseshoe is destroyed by a type one bifurcation that is homoclinic in the 
orientation preserving case, and heteroclinic otherwise. In either case we 
conjecture that this bifurcation is the first bifurcation among all orbits of 
the Henon map as we recede from the anti-integrable limit. 

With our continuation method, we are able to assign a "global code" to 
each orbit, by fixing the designation to that at the AI limit. In the Henon 
map, we demonstrate that this AI code is equivalent to the standard horse- 
shoe code (when it exists), but it also gives a consistent way of assigning 
symbols to orbits beyond the destruction of the complete horseshoe. Re- 
markably, there appears to be a relationship between the AI codes for a 
number of systems including billiards, twist maps of the cylinder and the 



Henon map. We will explore this relationship in a forthcoming paper |29|. 

We relate the properties transition time, type, and signature of homo- 
clinic orbits to properties of the core sequence. We also demonstrate that the 
ordering of the homoclinic orbits on the manifold segments U and S is the 
standard unimodal ordering. The notion of double neighbors and lemmas ^ 
and P give a necessary condition for a pair of homoclinic orbits to bifurcate. 
Surprisingly, these also give a forcing relation that tells us which homoclinic 
bifurcations have to occur before other ones. Showing that homoclinic bi- 
furcations can only take place between double neighbors with the same core 
length. Lemma |ll| gives a symbolic criterion for a pair of homoclinic orbits 
to be neighbors. The ordering is certainly valid until the complete horseshoe 
is destroyed, which leads to the theorem that the first homoclinic bifurca- 
tion of the hyperbolic fixed point in the area preserving Henon map occurs 
between the pair of type one orbits. When b = ±1, the Henon map has a 
symmetry and we discuss the mechanism by which pitchfork and asymmet- 



are 



ric saddle node bifurcations occur. The key ingredients to Theorem 16 
the ordering on the manifolds and the existence of a reversor for the map. 
As a result, the scenario of a tip of the manifold just repeatedly piercing 
through the other manifold (which is most natural in the case without sym- 
metry) is impossible. Among the possible alternatives are the occurrence of 
a pitchfork bifurcation or the creation of an asymmetric saddle-node pair, 
which are both not generic in the non-reversible case. 

With our continuation technique we compute numerical values for var- 
ious bifurcations of the homoclinic orbits up to type eleven. We sketch 
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the bifurcation diagram at type three, and then use a simple algorithm to 
construct the much more complex figures for higher core length. 

In contrast to our method for finding periodic orbits, the Biham and 
Wenzel method |24, ^] is known to fail in certain cases and can only be 
justified in the neighborhood of the AI limit [11]. Nevertheless, for the area 
preserving Henon map, we observe precisely the same number of orbits in 
the main DMS gap using our technique as was reported by DMS using the 
Biham- Wenzel method . We extend the original experiment of DMS in two 
ways. First, we study an order of magnitude more orbits than the original 
experiment and yet the gaps originally reported by DMS persist. Second, we 
observe that homoclinic bifurcations are responsible for these gaps and we 
list the symbolic labels of the orbits that form the gap endpoints in Table 
p. These gaps correspond to the creation and destruction of parameter 
intervals where the dynamics of the area preserving Henon map appears to 
be conjugate to a subshift of finite type. We find a similar gap structure 
for a particular subshift, Sjf, and list the symbolic labels for the homoclinic 
orbits that form the endpoints of the gap that is the analog of the DMS 
gap. The role of the special symmetric saddle- node bifurcations without an 
accompanying pitchfork bifurcation in this scenario remains to be elucidated. 

Many of our analytical results could be transferred from the area pre- 
serving case 6 = 1 to the orientation reversing case b = —1. In this case 
there also exists a reversor, so that the theorem about the first bifurcation 
and about impossibility of pitchfork bifurcations could be generalized to this 
case. The main difference is that now we are not studying homoclinic orbits 
of a fixed point which is invariant under the reversor, but instead hetero- 
clinic orbits connecting fixed points that are mapped into each other by the 
reversor. Correspondingly the action of the reversor on symbol sequences is 
not just reading them backwards, but it is reading backwards and flipping 
each symbol. 

Let us finally remark that we can extend the antimonotonicity result 
[13| to the non-dissipative case b = ±1 because there are bifurcations that 
do occur in the "wrong" direction, i.e., orbits that are created when k is 
decreased. This is described in detail in a separate paper [0] where it is 
related to the fact that in the neighborhood of the period tripling bifurcation 
the twist generically vanishes. 
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